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Development  of  a Corona  Anemometer 
for  Measurement  of  Stratospheric  Turbulence 


1.  INTRODUCTION 

Man  has  the  potential  for  discharging  pollutants  into  the  stratosphere  that  alters 
its  chemical  balance  and  results  in  serious  environmental  harm  to  the  earth's  cli- 
mate and  living  organisms.  This  threat  has  been  recognized  and  national  efforts 
are  underway  to  determine  what  harmful  effects  can  be  generated  by  the  release  of 
chlorofluoromethanes,  principally  from  aerosols,  release  of  NOx  in  aircraft  and 
missile  exhausts  and  other  sources  of  NOx  and  halocarbons.  The  steps  leading  to 
an  assessment  of  these  potentially  harmful  sources  involve  careful  examination  of 
stratospheric  composition  and  development  of  simulation  models  of  the  chemical 
and  transport  processes. 

Turbulence  is  an  important  parameter  describing  the  mixing  of  gases  in  the 
stratosphere.  How  rapidly  a polluting  gas  mixes  and  is  transported  vertically  indi- 
cates how  long  the  pollution  remains  in  the  stratosphere.  This  residence  time  is 
key  to  predicting  the  harmful  effect  a particular  pollutant  will  have.  Estimates  of 
one  to  two  years  residence  time  is  one  reason  why  very  small  concentration  levels 
of  pollution  can  produce  significant  changes  in  the  stratospheric  composition.  At 
the  present,  no  detailed  picture  exists  that  represents  the  actual  physical  processes, 
including  turbulence,  circulation,  and  winds,  existing  in  the  stratosphere.  This  has 
been  due  to  lack  of  realization  in  earlier  years  that  the  stratospheric  environment 
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could  be  altered  and  result  in  serious  changes  in  man's  world.  Until  recently  there 
have  been  too  few  measurements  conducted  to  determine  the  transport  processes. 
One  inhibition  to  measurements  has  been  the  absence  of  satisfactory  measurement 
techniques  for  stratospheric  turbulence.  This  report  describes  the  development 
and  testing  of  a corona  anemometer  adapted  to  measure  the  stratospheric  turbulence 
processes. 

2.  DEVELOPMENT  OF  CORONA  SENSOR 

The  corona  anemometer  for  detection  of  turbulence  (CADT)  measures  wind 
velocity  and  turbulence  by  sensing  the  displacement  fluctuation  of  ion  current 
on  a target.  The  CADT  sensor  consists  of  a corona  source  generating  a flow  of 
ions  from  the  corona  plate  to  the  target  plate.  This  is  illustrated  for  one-dimen- 
sion in  Figure  1.  The  ions  are  generated  in  the  electric  field  created  between  the 
high  voltage  applied  to  a fine  point  needle  and  the  voltage  applied  to  the  corona  con- 
trol plate.  A small  hole  in  the  control  plate,  1/4  fn.  diameter,  on  the  centerline  of 
the  needle  allows  a fraction  of  the  corona  current  «10  percent  to  be  drawn  to  the 
target  plate.  The  target  plate  is  segmented  into  halves  at  each  side  of  the  center- 
line.  The  ion  stream  diffuses  outward  as  it  transverses  the  air  gap  between  the 
two  plates.  In  the  absence  of  a wind,  the  current  collected  by  each  half  of  the  tar- 
get plate  is  identical.  With  a wind  blowing  between  the  two  plates,  the  ions  are 
deflected  downstream  by  the  wind.  The  difference  in  current  collected  on  the  two 
plates  is  a measure  of  the  wind  speed. 


Figure  1.  Stratospheric  Corona  Anemometer  Schematic 
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The  advantage  of  CADT  for  stratospheric  wind  measurement  is  the  ease  by 
which  the  sensor  response  to  windspeed  can  be  maintained  uniform,  while  the  force 
of  the  wind  dramatically  changes  with  altitude  due  to  the  varying  atmospheric  den- 
sity. Mechanical  anemometers,  in  contrast,  are  optimized  to  a particular  density 
altitude. 

Measurement  of  stratospheric  turbulence  is  very  difficult.  The  stratosphere, 
in  contrast  to  the  troposphere,  has  a low  density  and  consequently  low  levels  of 
turbulent  energy.  This  means  that  the  usual  turbulence  detectors  employed  in  the 
troposphere  are  inadequate.  Mechanical  anemometers  have  too  much  inertia,  the 
energy  threshold  of  acoustic  detectors  is  too  high,  and  instrumental  errors  domi- 
nate the  output  of  hot  wire  detectors.  Optical  anemometers  require  a large  density 
of  light  scattering  particles  which  are  not  found  in  sufficient  density  in  the  strato- 
sphere. A new  type  of  turbulence  detector  was  needed  for  the  stratosphere. 

A suitable  anemometer  is  one  based  on  measuring  the  advective  motion  of  ions 
traversing  between  two  opposed  plates.  The  ions  can  be  produced  by  a stable  or  co- 
rona with  electrical  control  of  the  corona  current  to  maintain  a uniform  flow  of  ions 

The  corona  anemometer  has  been  used  to  study  turbulence  in  wind  tunnels*  and  was 

2 

first  applied  for  stratospheric  turbulence  measurement  by  Barat.  A corona  ane- 

3 

mometer  to  measure  gas  flow  rates  was  successfully  developed  by  Liliet.i'ield  for 
use  on  stratospheric  particle  samples.  This  report  describes  the  development  of 
a corona  anemometer,  its  calibration,  and  the  results  of  the  initial  stratospheric 
turbulence  measurements. 


2.1  Theory  of  Operation 

A corona  is  an  electrical  breakdown  of  a gas  which  occurs  where  nonuniform 
electric  fields  are  encountered.  It  can  occur  with  asymmetrical  electrode  geome- 
tries such  as  point  to  plane,  sphere  to  cylinder,  or  wherever  sharp  edged  or  small 
radii  electrodes  appear.  Breakdown  is  initiated  by  the  high  fields  concentrated  at 
one  or  both  of  the  electrodes  but  complete  breakdown  to  an  arc,  at  or  near  the 
corona  threshold,  is  precluded  by  the  weak  field  between  the  electrodes.  The  ap- 
pearance of  corona  current  at  the  threshold  levels  is  usually  unstable  due  to  an 
insufficient  number  of  triggering  electrons.  Above  threshold  and  at  a current  on 
the  order  of  microamperes,  steady  corona  is  achieved.  Here  an  increase  in  curren 


1.  Franzen,  B.,  Fucks,  W. , and  Schmitz,  G.  ( 1961)  jCorona  anemometer  for 

measurement  of  turbulence  components,  Zelt.  fur  Flugurssenschaften 
9:347-351. 

2.  Barat,  J.  (1975)  Une  methode  de  mesure  directe  du  taux  de  disspation  d'energie 

turbulente  dans  la  stratosphere.  Compt.  Rend.  Acad.  Sci.  Paris  B281: 
53-56. 

3.  Lilienfeld,  Pedro,  Solon,  L. , and  DiGlovanni,  Hugo  (1967)  Ion  tracer  anemom- 

eter for  the  measurement  of  low  density  flows.  Rev.  Sci.  Inst.  38:405-409. 


is  nearly  proportional  to  a potential  increase.  This  is  the  so-called  Ohm's  Law 
Regime  in  which  the  CADT  sensors  operate  most  frequently.  Beyond  this  region, 
current  increases  parabolically  with  potential  and  eventually  results  in  a complete 
transient  arc  or  spark  breakdown.  Tests  performed  with  CADT  sensors  reveal 
that  steady  corona  onset  is  achieved  at  tens  of  volts  above  threshold  and  spark 
breakdown  occurs  at  a few  hundred  volts  above  onset.  Actual  values  depend  upon 
pressure  P,  and  electrode  separation,  / . Threshold  voltages  follow  the  well  known 
Paschen  Law  where  Vg  = Vg  (P /). 

The  basic  theory  of  operation  of  the  corona  anemometer  can  be  predicted  from 
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the  corona  theories  developed  by  Loeb,  * Brown,  and  McDaniel  and  Mason. 

Corona  anemometers  have  been  developed  to  measure  turbulence  in  wind  tunnels1 

and  to  measure  surface  winds.  A new  development  was  required  to  adapt  the 

corona  sensor  to  measure  stratospheric  winds  relative  to  the  motion  of  a floating 

balloon  and  to  directly  detect  the  existence  of  turbulence  and  measure  the  turbulent 

intensity. 

Conceptually,  the  displacement  of  the  ions  can  be  seen  as  the  result  of  the 
interaction  of  the  ion  stream  and  airstream.  The  displacement  angle,  0,  is  defined 
by  the  normal  displacement  length  D,  of  the  ion  stream  from  the  centerline  and  by 
the  path  length  L over  which  the  ions  drift  to  the  target.  For  small  deflections, 
conservation  of  linear  momentum  requires  that  the  tangent  of  fl  be  equal  to  the  ratio 
of  the  initial  momenta  of  the  airstream  to  ion  stream.  This  is  represented  as: 


tan  0 


D 

L 


(mU)  . 
air 

^Vions 


(1) 


where  m is  the  mass  of  air  or  ions  which  interact  within  the  ion  stream  volume  V; 
U is  the  freestream  air  velocity,  and  v^  is  the  ion  drift  velocity. 

To  a first  approximation  the  momentum,  mil,  can  be  related  to  the  mass  flow 
(pUA)by  accounting  for  the  interaction  volume  V and  the  respective  cross  sectional 
area  of  interaction,  A,  where: 


4.  Loeb,  L.  B.  (1960)  Basic  Processes  of  Gaseous  Electronics.  Univ.  of 

California  Press,  Los  Angeles. 

5.  Loeb,  L.  B.  (1965)  Electrical  Coronas.  Univ.  of  California  Press,  Los 

Angeles. 

6.  Brown,  S.C.  (1966)  Introduction  to  Electrical  Discharges  in  Gases.  John 

Wiley  & Sons,  New  York. 

7.  McDaniel,  E.W.,  and  Mason,  E.A.  (1973)  The  Mobility  and  Diffusion  of 

Ions  in  Gases,  John  Wiley  & Sons,  New  York. 

8.  Waletzko,  J.A.  (1975)  A new  ion  displacement  system  to  measure  the  two 

dimensional  wind  vector,  Tech.  Conf.  on  Automated  Meteorological  Systems. 
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mU  = (pUA)  x J 


(2) 


V 


It  follows  that, 

JD.  (pUA)airXV/Aair 
L (PivdA)ions  * V/AT 


(3) 


where  AT  is  the  area  of  the  target  upon  which  the  ions  impinge.  Experiments  have 
shown  that  by  maintaining  the  potential  difference  Ep  constant  between  control  plate 
and  target,  AT  is  a weak  function  of  p.  Also  the  total  current  collected,  IT>  is  pro- 
portional to  the  ion  mass  flow,  IT  oc  (pvdAT)ions.  This  results  in  the  ability  to 
measure  mass  flux  of  the  airstream  as: 


p u = C 1 D IT 

where  Cj  is  a constant  containing  the  physical  parameters  of  the  sensor. 

The  details  of  operation  in  terms  of  ion  mobility  and  corona  current  voltage 
relationships,  have  been  described  by  Waletzko.  The  corona  ions  acquire  a drift 
velocity  parallel  to  the  electric  field  across  the  two  plates  as 

V(J  = n • E (4) 

where  p is  the  ion  mobility  and  E is  the  electric  field  established  by  the  potential 
difference  between  the  two  plates  separated  a distance  L.  This  linear  relation  holds 
as  long  as  E/p  < 2 volt/cm  torr  and  if  E is  uniform.  The  ion  transit  time  between 
the  two  plates  is  therefore 


L 


o 


(5) 


The  distance  the  ions  are  swept  downstream  during  the  transit  time  is 


D = U . r 
air 


dZ 

ME  • 


(6) 
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The  ion  mass  and  the  air  molecule  mass  are  nearly  identical  so  that  the  ions 
acquire  a velocity  component  parallel  to  the  wind  due  to  momentum  transfer.  This 
velocity  component  is  acquired  in  a time  period  of  a few  collisions  (t  = 1-10  nsec) 
which  is  small  in  comparison  to  the  ion  transit  time. 

The  corona  current  is  collected  on  an  electrode  consisting  of  a matrix  of  small 
resistors.  The  resistors  are  connected  together  in  such  a manner  as  to  effectively 
measure  the  current  falling  on  each  of  segment.  Subsequent  development  and  flight 
experience  has  shown  the  multi-element  matrix  to  have  insufficient  velocity  resolu- 
tion due  to  the  spacing  between  elements.  A four  plate  electrode  is  currently  being 
used  with  excellent  results.  Considering  only  one -dimens  ion,  the  tisplacement  of 
the  ion  beam  is  sensed  by  measuring  the  difference  in  current  collected  on  the 
plates  (see  Figure  2).  The  corona  current  collected  has  been  measured  across  the 
surface  of  a test  electrode  and  found  to  have  a Gaussian-shaped  distribution  with  a 
half  peak  width  of  =5  mm.  With  this  symmetric  distribution,  the  difference  in  cur- 
rents represents  a displacement  of  the  ion  beam  and  can  be  visualized  as  a differ- 
ence in  beam  area  projected  on  each  half  of  the  electrode.  For  a uniform  current 
distribution  the  area  difference  in  terms  of  displacement,  D,  is: 


(7) 


where  r is  the  radius  of  the  corona  beam.  This  can  be  approximated  using  the  fol- 
lowing expression  as  shown  in  Figure  3a 


Figure  2.  Corona  Target  Plates 
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For  a narrow  corona  ion  beam  of  radius  r,  the  difference  in  the  currents  is  pro- 
portional to  the  displacement  of  the  corona  beam  from  the  electrode  center. 


Aq  " A i 

l2  ‘ *1  3 A *T 


(9) 


For  small  displacements,  the  current  difference  is  linearly  proportional  to  the  dis- 
placement 


h-h-U?)'  T 


(10) 


The  corona  beam  radius  is  subject  to  change  due  to  both  lateral  diffusion  and  effects 
of  the  nonparallel  electric  field,  that  is.  r = L / <Ep>.  If  only  laterial  diffusion 
existed,  the  beam  radius  function  would  be  (McDaniel  and  Mason,  1973) 


/j(  Ep)  0.173(E„L) 


-1/2 


(11) 


The  expression  relating  wind  velocity  to  a measurement  of  corona  current  can  now 
be  written.  Combing  Eqs.  (6)  and  (10),  the  velocity-current  expression  is 


I2  * It 


7 LJ"V  ' --■> 


/ dz 
J U E. 


U 


(12) 


The  ambient  mobility  of  the  positive  ion,  n.  >9  defined  in  terms  of  the  reduced 


mobility,  Kq  and  gas  density  as 


(13) 


The  reduced  mobility  is  observed  to  be  independent  of  the  electrical  field  provided 

7 

E/p  < 30  (v/cm  torr)  (McDaniel  and  Mason  ). 
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Thus  the  integral  in  Eq.  (12)  can  be  evaluated  as 
L 


(14) 


with  the  result  of  current  difference  given  as 


i t - 4 g<Z>  o ^ Uair 
2 1 ' *o  WP  ' 


(15) 


When  the  CADT  was  operated  with  fixed  target  current,  l_,  the  current  difference 

/2 


was  observed  to  be  a function  of  p U . and  having  an  E V*  dependence.  This  im- 

r 1/2  alr  P 

plies  that  i (E  ) = /3  E ' and  3 is  a proportionality  constant  as  shown  in  Eq.  (11) 


J 


'T’U“ 


air 


(16a) 


The  CADT  can  be  operated  at  a fixed  E^  field  strength  with  a mass  flow  detection 


(16b) 


»K  p ed/<ed) 

where  the  constant  C = 4^g(z) *3  determined  by  calibration.  This  is  the 


basic  velocity  relation  used  to  interpret  measurements  with  the  corona  anemometer 


operated  at  fixed  1^.  and  Ep. 


2.  1.  1 ION  GENERATION 


The  corona  can  be  generated  in  many  ways:  sphere -to-cylinder,  point-to-point 
or  point -to-plate.  In  this  instrument,  the  corona  ions  are  created  by  high  voltage 
breakdown  between  a fine  needle  point  separated  3 mm  from  the  hole  in  the  upper 
plate.  The  gap  distance  between  the  needle  tip  and  the  hole  is  adjusted  until  a fixed, 
2 pA  current  exists  when  1700  volts  is  applied  at  120  mm  Hg  pressure  in  a dry  air 
atmosphere.  The  corona  current  generated  can  be  expressed  as 


Ic  = (V  - Vg>  vd  h (17) 

where  V is  the  potential  difference  between  thtf  needle  and  plate,  Vg  the  threshold 
voltage  for  creation  of  corona  current  (Geige.  counter  region)  and  h is  a complex 


geometric  parameter  defining  the  two  electrodes.  At  the  high  field  strength 
E/p  > 2 V/cm-torr,  the  drift  velocity  follows  the  relation 


v . 

d o 


The  corona  current  relation  when  the  applied  voltage  V is  much  greater  than 


I = (K  h)(V 
c o 


where  1 represents  the  separation  of  approximately  3 mm  between  the  needle 
point  and  the  plate.  The  corona  instrument  is  operated  to  maintain  a constant 
corona  current  throughout  the  flight.  Thus  it  is  necessary  to  lower  the  high  volt- 
age potential  on  the  needle  as  the  instrument  ascends  to  higher  altitudes.  Labora- 

' 1/2 

tory  measurements  at  constant  temperature,  shown  in  Figure  4,  confirm  the  p ' 
dependence  of  the  corona  current. 

A small  fraction  of  the  corona  current  generated  is  drawn  through  the  6 mm 
hole  in  the  upper  plate,  down  to  the  collector  plate.  The  fraction  of  current 


< 1000 


(V-VgJv'V/Ic 
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Figure  4.  Corona  High  Voltage  Dependence  on  Pressure 
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drawn  (=10  percent)  depends  upon  the  ion  density  in  the  cavity  region  and  the  field 
strength  through  the  hole.  The  total  current  collected  is  referred  to  as  the  target 
current  Lj,.  The  current  collected  has  a nonlinear  dependence  on  E^  and  a weak 
linear  dependence  on  1^.  The  total  target  current  increases  with  pressure  to  a 
pressure  between  5 and  20  mm  (depending  upon  E^),  after  which  the  target  current 
decreases  and  eventually  reaches  zero. 

Detailed  observations  reveal  that  the  target  current  is  independent  of  wind 

-5  3 

velocity  at  air  densities  below  2X10  g/cm  . At  higher  densities,  the  target 
current  decreases  with  increasing  velocity.  This  decrease  in  current  is  believed 
to  occur  in  the  neighborhood  of  the  upper  plate  hole  where  a portion  of  the  current 
is  blown  to  the  upper  plate  by  air  passing  up  into  the  cavity. 

2.2  Engineering  Design 

Engineering  development  of  the  sensor  components  required  evaluation  and 
optimization  of  both  the  electrical  performance  and  aerodynamic  performance. 

Wind  tunnel  testing  was  used  to  evaluate  the  aerodynamic  performance  (see  section 
2.  4)  and  select  an  optimized  geometry.  The  first  parameter  to  be  chosen  was  the 
plate  separation  distance  L.  A 2.  5 cm  (1")  separation  was  chosen  as  a compromise 
between  sensitivity,  smooth  airflow  and  minimization  of  spread  of  the  corona  ion 
beam  across  the  gap.  A large  radius  of  curvature  at  the  edges  of  the  plate  was 
found  to  be  necessary  to  prevent  separation  of  flow  at  the  edges  which  would  lead 
to  internally  generated  turbulence.  The  stand-off  rods  holding  the  top  plate  to  the 
sensor  were  found  to  be  useful  as  a means  of  calibrating  the  flow.  The  wake  behind 
the  rods,  under  certain  Reynolds  number  conditions,  has  regularly  spaced  eddies 
which  produce  a single  frequency  as  they  are  carried  across  the  corona  beam  at  the 
speed  of  the  airstream.  The  airstream  velocity  can  be  measured  by  observing  the 
frequency  of  arrival  of  the  turbulent  eddies. 

Fabrication  and  shaping  of  the  high  voltage  probes  was  a key  factor  in  the  suc- 
cessful development  of  CADT.  The  high  voltage  probe  must  have  a sharp  point  to 
achieve  steady  operation  of  a corona  at  voltages  considerably  below  spark  break- 
down voltages  and  to  minimize  the  probe  heating.  Heat  conduction  from  the  point 
itself  to  a heat  sink  is  vital  in  order  to  preserve  the  integrity  of  the  sharpened  probe 
and  achieve  20-30  hour  operational  lifetimes.  The  composition  of  the  needle,  there- 
fore, must  be  such  that  the  probe  has  rigidity,  is  machinable  and  has  good  conduc- 
tivity. Numerous  materials  were  tried  and  rejected  before  finding  the  best  mate- 
rial. Tooled  steel  (sewing  machine  needles)  had  poor  conductivity  and  the  points 
rapidly  melted.  Tungsten  carbide  was  too  hard  to  form  an  adequate  point.  Gold 
and  platinum  did  not  last  long.  A needle  formed  from  90  percent  platinum  and 
10  percent  Iridium  was  found  to  hold  a good  point  in  excess  of  20  hr  of  operation. 
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The  best  design  was  found  using  a 3 mm  long,  0.  6 mm  diameter  needle  with  a 
20  degree  taper  to  a tip  diameter  of  0.  08  mm.  The  needle  was  found  to  survive 
the  occasional  arc-overs  that  occurred  during  various  tests.  An  excessive  number 
of  hours  of  operation,  however,  results  in  a ball  forming  on  the  tip  of  the  needle. 

The  high  voltage  probe  is  mounted  in  a tapered  holder.  The  taper  =60  degree 
included  angle  is  required  to  achieve  the  best  heat  conduction  by  getting  the  copper 
holder  as  close  to  the  probe  point  without  producing  arc -over.  The  needle  is  initially 
polished  to  a fine  tapered  (=20  degree  angle)  point.  The  CADT  is  operated  for  two 
hours  to  achieve  a "burn  in"  of  the  probe  point.  Initially,  the  high  voltage  required 
for  corona  is  about  10  percent  lower  than  the  final  operating  level.  The  probe  point 
during  the  first  hours  of  operation  undergoes  melting,  and  melts  back  to  a small 
radius  point  =0.  1 mm,  requiring  a slightly  higher  voltage  for  corona.  Thereafter, 
the  voltage  remains  constant.  The  anticipated  velocities  range  between  0 and  20 
m/sec  with  a resolution  requirement  of  1 cm/sec.  The  velocity  range  that  can  be 
measured  is  a function  of  the  ion  drift  velocity.  This  is  determined  by  the  choice 
of  the  plate  voltage,  Ep.  A high  Ep  implies  a fast  drift  velocity  and  a small  deflec- 
tion for  a given  airstream  velocity.  A constant  Ep  of  350  volts  was  found  to  give 
the  desired  resolution  and  operate  the  corona  in  the  desired  linear  range  over  the 
expected  velocity  range.  (The  plate  voltage  is  correctly  Ep  • L.  ) 

The  CADT  operates  within  two  constraints.  Too  high  an  operation  potential 
creates  a strong  field  between  the  needle  and  top  plate,  reduces  the  current  1^. 
that  can  be  pulled  to  the  collector  plate.  Decreasing  the  applied  potential  in- 
creases the  collected  current  until  the  applied  potential  approaches  the  corona 
threshold  potential,  a point  where  no  corona  will  occur.  Thus  operating  the  high 
voltage  (HV)  to  produce  the  fixed  2 pA  of  corona  current,  requires  dropping  the 
Ep  at  pressures  below  15  mm  Hg  to  100  volts  in  order  to  be  above  the  threshold 
potential.  This  increases  the  velocity  sensitivity,  that  is,  a given  velocity  requires 
a larger  AI/I^.  for  a smaller  Ep.  This  limits  the  maximum  velocity  that  can  be 
detected  before  Al/1^.  approaches  unity.  Subsequent  flight  experience  has  shown 
the  velocities  to  be  lower  at  the  high  altitude  so  that  dropping  Ep  is  not  expected 
to  cause  problems. 

The  high  corona  voltage,  Ecf,  required  to  produce  the  desired  constant  corona 
current,  I varies  with  the  atmospheric  density  as  shown  by  Eq.  (19)  and  measured 
in  Figure  4.  A circuit  was  designed  to  vary  the  high  voltage  as  required  to  achieve 
a constant  corona  current  of  two  pA  with  a 12  sec  response  time.  The  high  voltage 
ranges  from  1700  volts  at  120  mm  Hg  to  700  volts  at  15  mm  Hg  pressure.  The 
electric  field  between  the  two  plates  is  perturbed  by  the  high  voltage  which 
creates  a small  nonlinear  variation  in  velocity  response  with  altitude.  Calibration 
of  the  CADT  is  required  over  the  anticipated  operating  region.  Heads  have  been 
assembled  from  identical  components  in  a very  precise  manner.  Variations  in 
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performance  between  "identical"  units  usually  can  be  traced  to  some  slight  mis- 
alignment that  can  be  corrected. 

The  output  of  the  CADT  sensor  is  measured  on  each  axis.  Four  data  channels 
are  required  to  measure  the  high  and  low  gain  outputs  of  the  x and  y axis  with  the 
high  gain  output  ten  times  the  low  gain.  These  channels  measure  primarily  the 
mean  windspeed.  The  turbulence  is  measured  with  four  data  channels  having  a 
high  pass  filter  removing  signals  below  20  Hz.  There  are  four  ac  channels,  ampli- 
fied a factor  ten  over  the  dc  channels,  which  measure  the  velocity  fluctuations  be- 
tween 20  Hz  and  2000  Hz. 

The  developed  sensor  is  shown  in  Figure  5.  A total  of  four  sensors  were  built, 
with  the  best  two  flown  in  the  initial  flight  test,  described  in  Section  3. 


Figure  5.  Flight  Ready  Photograph  of  Corona  Anemometer 

2.  2.  1 ELECTRONICS  DESIGN 

The  current  from  each  target  A,  B,  C,  and  D goes  to  the  preamp  (see  Figure  6) 
where  it  is  converted  to  a voltage,  through  the  preamp,  and  to  the  output,  where  the 
output  is  low  impedance.  The  preamp  also  sums  all  of  the  current  that  is  collected 
on  the  four  targets  and  this  is  called  E,^  and  is  measured  as  a voltage  across  a 
resistor.  To  measure  a signal  in  the  +x  direction,  the  x amplifier  adds  the  volt- 
ages as  A + B and  C + D,  then  subtracts  the  sum  of  (C  + D)  from  (A  + B).  This 
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latter  signal  (Ey)  is  then  sent  to  a function  card  where  it  is  divided  by  the  product 
of  and  a voltage  proportional  to  the  ambient  pressure.  The  pressure  correc- 
tion circuit  corrects  for  the  offset  and  scale  factor  of  the  pressure  transducer  and 
produces  a linear  voltage  output.  The  signal  Exj(Ej • p)  is  then  sent  back  to  the 
amplifier  where  various  gained  outputs  are  produced.  The  dc  output  is  actually 
the  signal  with  only  a high  frequency  cut  off,  the  ac  outputs  goes  through  an  addi- 
tional low  pass  filter  leaving  10  Hz  to  1 KHz  with  a 3 db  roll  off  on  the  high  end. 

From  the  y amplifier,  the  output  is  fed  to  an  RMS  module  and  the  RMS  value 
is  then  transmitted. 

A zero  correction  circuit  is  used  to  correct  for  two  errors,  one  being  the 
long  time  drift  of  the  amplifiers  and  the  more  important  reason  being  the  change 
in  beam  pattern  with  time  due  to  wear  of  the  point.  When  the  shutter  is  closed,  a 
cylindrical  shell  encloses  the  sensing  head  air  gap.  A check  is  made  of  the  output 
of  the  dc  low  gain  amplifier.  If  it  is  not  at  zero,  a digital  counter  generates  a volt- 
age output  until  the  output  amplifier  is  zero.  At  that  point,  it  stops  counting  and 
holds  this  offset  until  it  is  called  for  again. 

The  corona  current,  the  current  from  point  to  control  plate,  is  maintained  at 
2 ^A  ± 0.03  for  all  pressures.  The  corona  current  goes  through  an  isolation 
amplifier,  the  output  of  the  amplifier  is  a linear  function  of  current  and  this  is  fed 
to  an  electronic  servo  High  Voltage  Regulator.  The  output  is  fed  to  a dc  to  dc 
high  voltage  power  supply  and  maintains  constant  current.  The  high  voltage  con- 
trol is  the  voltage  from  the  high  voltage  regulator  and  this  is  a function  of  pressure. 
This  control  voltage  is  used  to  switch  Ep  from  350V  at  relatively  high  pressures  to 
100V  at  low  pressures.  When  we  shift  the  voltage,  we  also  change  the  gain  of 
the  E,^  gain  control.  The  commutator  is  a multiplexer  used  to  sample  all  of  the 
different  functions  necessary  to  determine  if  the  system  is  operational.  These 
sampled  voltages  would  be  the  output  of  B+  and  B-  regulators,  battery  voltages, 
shutter  position  open  or  closed,  main  boom  position,  and  other  housekeeping  data. 
Figure  7 shows  the  electronic  block  drawing  in  greater  detail. 

2.3  Calibration 

The  CADT  system  was  calibrated  for  velocity  using  two  approaches.  The  first 
calibration  was  conducted  inside  a 45  cm  (18  in. ) laboratory  bell  jar  using  a small 
turbine  air  blower  to  generate  the  velocity  through  the  sensor.  The  second  calibra- 
tion was  performed  at  the  Wright  Brothers  Low  Speed  Wind  Tunnel  at  MIT. 

The  bell  jar  calibration  was  set  up  as  shown  in  Figure  8.  The  flow  velocity 
inside  the  duct  was  measured  using  a Pitot  tube  and  Baratron  absolute  and  differ- 
ential pressure  sensors.  The  velocity  calibration  is  given  in  Figure  9.  The  CADT 

system  velocity  calibration  in  terms  of  voltage  output  V or  V is  presented  in 

x y 
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Figure  10.  Separate  calibrations  for  each  sensor  and  the  X and  Y directions  are 
required  but  they  are  similar  in  response.  The  velocity  in  each  axis,  U and  U 

x y 

is  determined  to  be  functions  of  the  measured  voltage  in  each  axis. 


U = V 
x x 


U = V 

y y 


Low  gain 


ux  = Vx/10 


U = V /10 

y y 


High  gain 


The  velocity  calibration  could  not  be  accurately  conducted  at  speeds  less  than 
0.  5 m/sec  at  low  pressures  due  to  the  lack  of  an  accurate  and  independent  velocity 
standard.  A test  of  the  linearity  of  the  velocity  response  was  developed  and  it  demon- 
strated the  CADT  has  excellent  linear  velocity  response.  Thus  the  calibration  con- 
sists of  performing  an  accurate  measurement  over  one  velocity  range  a. id  establish- 
ing the  existence  of  a linear  response  over  the  entire  velocity  range.  The  linearity 
was  proven  using  a special  plate  that  permitted  translating  the  collector  plate  along 
the  velocity  vector  in  steps  as  small  as  0.  001  in.  The  procedure  was  to  establish 
a given  velocity  and  measure  the  deflection,  AI/Irp.  The  flow  was  turned  off  and 
the  collector  plate  translated  relative  to  the  needle  until  the  same  AI  1^.  was  ob- 
served. The  flow  is  then  turned  on  and  the  deflection  2 AI/I,j.  is  measured.  The 
flow  is  then  turned  off  and  the  collector  plate  translated.  This  was  repeated  in 
increments  out  to  a AI/I  0.33  as  shown  in  Figure  3b. 

A calibration  in  flight  can  be  obtained  by  taking  advantage  of  the  existence  of 
vortex  shedding  from  the  standoffs  separating  the  two  plates.  A regular  spaced 
Karman  vortex  street  is  generated  as  the  uniform  flow  separates  passing  around 
the  circular  cylinder  standoffs  separating  the  two  plates.  The  shedding  frequency, 
n,  of  the  individual  vortexes  can  be  made  dimensionless  in  terms  of  the  cylinder 
diameter,  D,  and  the  flow  velocity,  U,  as 


where  S is  shown  as  the  Strouhal  number.  The  vortex  is  observed  to  depend  on  the 
flow  Reynolds  number.  The  vortex  street  is  stable  and  regular  only  between 
Reynolds  No.  40  and  150  (Roshko,  ''  Webster,  10  and  Tritton,  ll).  The  vortex  street 


9.  Roshko,  R.  (1954)  On  the  Development  of  Turbulent  Wi.kes  From  Vortex 
Streets,  NACA  Rep  1191. 

10.  Webster,  C.A.G.  (1954)  An  experimental  study  of  turbulence  in  a density- 

stratified  shear  flow,  Fluid  Mech.  19:221-245. 

11.  Tritton,  D.J.  (1970)  F.xperiments  on  the  Flow  Past  a Circular  Cylinder  at 

Low  Reynolds  Numbers. 
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Figure  8.  Bell  Jar  Velocity  Calibration  System 
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Figure  9.  Air  Velocity  as  a Function  of  Bell  Jar  Blower  RPM  and 
Air  Pressure 
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Figure  10.  Anemometer  Voltage  Calibration  as  a Function  of 
Airstream  Velocitv 


is  not  observed  at  Reynolds  number  below  40  and  the  Sfrouhal  number  is  a constant 
0.212  at  Reynolds  number,  Re,  greater  than  1000.  Below  Rp  150,  the  frequency 
occurs  from  instabilities  that  build  up  in  the  attached  laminar  wake.  At  Reynolds 
number  above  about  300,  the  frequency  occurs  from  instabilities  in  the  flow  separa- 
tion off  the  cylinder  itself. 

. 9 

The  Strouhal-Reynolds  number  relation  was  empirically  determined  by  Roshko") 

to  an  accuracy  of  1 percent  for  the  stable  region  as 


S = 0. 212 


For  the  turbulent  instability,  the  relation  is  found  as 


S = 0.  212  ^ 1 - j R{ 


The  transition  to  turbulent  flow  separation  occurs  earlier  in  the  presence  of  turbu- 
lent flow  when  the  velocity  fluctuations  exceed  0.  1 percent  (Webster  u). 

The  flow  velocity  as  a function  of  Strouhal  frequency  can  be  determined  from 
Eqs.  (21)  and  (22)  using  the  definition  of  Reynolds  number  and  expressed  as 


V 


(24) 


= _n2_  + 2L1 

“(cm/s)  0.212  D 

where  v is  the  kinematic  viscosity.  This  relationship  is  plotted  in  Figure  11  for 
various  values  of  viscosity. 

The  CADT  instrument  ac  signal  is  frequency  analyzed  to  evaluate  the  power 
spectrum.  A typical  spectrum  is  shown  in  Figure  12. 

The  flight  data  was  sampled  and  comparison  between  the  velocity  determined 
from  the  Strouhal  frequency  using  Eq.  (24)  and  the  velocity  determined  from  the 
corona  displacement  using  Eq.  (16).  The  results  are  shown  in  Figure  13.  Excel- 
lent agreement  is  found  between  the  two  velocities.  Note,  the  Strouhal  velocity  is 
resolved  into  x and  y coordinated  for  comparison  purposes  using  the  angle  of  the 
flow  determined  by  the  corona  sensor.  The  velocity  range  that  can  be  calibrated 
with  the  Strouhal  frequency  at  a given  altitude  is  shown  in  Figure  14. 

2.  3.  1 WIND  TUNNEL  INVESTIGATION  OF  CORONA  SENSOR 

An  instrument  to  be  placed  in  a moving  airstream  should  be  carefully  studied 
to  examine  what  perturbations  are  created  in  the  flow  by  the  instrument  itself. 

Ideally,  the  aerodynamic  investigations  should  be  conducted  at  simulated  tempera- 
ture, pressure  and  airspeed  conditions,  that  is,  full  scale  Reynolds  number.  The 
aerodynamic  tests  on  CADT  were  conducted  at  the  Wright  Brothers  Low  Speed 
Wind  tunnel  at  MIT.  This  tunnel  was  selected  for  its  capabilities  of  generating  low 

speed  and  low  turbulence  flows.  The  results  of  the  test  are  available  in  full  detail 
12 

(Durgin  and  Fanucci  ) with  a synopsis  presented  below. 

The  wind  tunnel  operates  only  at  atmospheric  pressure  and  temperature.  The 
stratospheric  corona  sensor  was  not  designed  to  operate  at  pressures  above  120  mm 
Hg  since  the  high  voltage  was  limited  to  2000  V.  Consequently,  it  was  necessary  to 
use  the  commercial  Thermal  Systems  Inc.  (TSI)  anemometer  in  the  wind  tunnel 
tests.  The  electronics  are  different  between  the  AFGL  and  TSI  instruments,  and 
the  separation  gap  was  a half  inch  for  TSI  and  one  inch  for  the  AFGL  instrument. 

Both  instruments  used  the  same  collector  and  anode  plates.  It  was  believed  that 
differences  were  not  important  for  an  investigation  of  the  aerodynamic  effects  and 
the  possibility  of  instrumental  induced  turbulence.  The  Wright  Brothers  Wind  tun- 
nel at  the  Massachusetts  Institute  of  Technology  was  operated  at  atmospheric  pres- 
sure at  speeds  from  0.  5 m/sec  to  14  m/sec  with  controlled  amounts  turbulence. 
Turbulent  flows  were  generated  by  installing  a wooden  grid  at  the  beginning  of  the 
test  section.  The  grid,  constructed  of  a square  lattice  of  1 in.  X 3 in.  boards  spaced 

12.  Durgin,  F.  H. , and  Fanucci,  J.  P.  (1977)  Static  and  Dynamic  Calibration  of  a 
Corona  Discharge  Anemometer,  AFGL-TR-77-0022. 
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Figure  11.  Strouhal  Velocity  Calibration 
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Figure  13.  Comparison  of  Bell  Jar  Velocity  Calibration  With  In- 
flight Calibration  Using  Strouhal  Velocity  Calibration 


Figure  14.  Altitude  Domain  Useful  for  Strouhal  Velocity  Calibration 

13 

12  in.  apart  (see  Figure  15)  was  designed  using  data  from  Baines  and  Peterson 
to  produce  a flow  which  varied  from  20  percent  turbulence,  2 m from  the  grid  to 
5 percent  turbulence  4 m from  the  grid.  Comparisons  of  spectra  taken  with  a hot 
wire  and  corona  anemometer  are  shown  in  Figure  16  through  19  for  a range  of  flow 
velocities.  It  is  evident  that  the  corona  anemometer  spectra  are  identical  to  hot 
wire  spectra  with  no  internally  generated  turbulence  present  at  these  turbulence 
levels.  At  very  low  levels  of  turbulence  generated  with  a special  fan,  turbulence 
is  generated  inside  the  corona  anemometer  as  shown  in  Figure  20.  It  is  apparent 
that  the  power  spectra  from  the  TSI  corona  anemometer  is  representative  of  free- 
stream  turbulence  when  the  freestream  turbulence  level  exceeds  one  percent.  This 
is  shown  in  Figure  21.  At  lower  turbulence  levels,  the  turbulence  generated  inside 
the  TSI  corona  anemometer  from  separated  flow  dominates  the  freestream  turbu- 
lence at  high  frequencies. 

The  CADT  must  accelerate  the  flow  passing  through  it  in  order  to  accommo- 
date the  boundary  layer  and  aperture  effect.  The  velocity  increase  was  measured 
with  the  hot  wire  inside  two  models,  each  having  a different  leading  edge  radius  of 
curvature.  Figure  22  indicates  the  flow  is  accelerated  10  percent  at  the  center  of 
corona  sensor.  The  amount  of  turbulent  fluctuations  inside  the  sensor  was  also 

13.  Baines,  W.  D. , and  Peterson,  E.  G.  (1951)  An  investigation  of  flow  through 
screens,  Trans.  ASME.  pp  467-480. 
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I Figure  lfi.  Power  Spectral  Density  Comparison  of  Hot  Wire  and  Corona 
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Figure  17.  Power  Spectral  Density  Comparison  of  Hot  Wire  and  Corona 
Anemometer,  U = 10,  6 m/sec,  X = 2m 
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Figure  18.  Power  Spectral  Density  Comparison  of  Hot  Wire  and  Corona 
Anemometer,  U = 2.  8 m/sec,  X 4,57 
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Figure  19.  Power  Spectral  Density  Comparison  of  Hot  Wire  and  Corona 
Anemometer,  U = 1.  6 m/sec,  X * 4.  57 
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Figure  20.  Power  Spectral  Density  Comparison  of  Hot  Wire  and  Corona 


VELOCITY  ACCELERATION,  U/U  AT  CENTER 

GO 


Figure  22.  Velocity  Acceleration  Produced  Inside  the  Anemometer 


measured  with  the  results  shown  in  Figure  23.  The  entrance  radius  of  curvature 
has  a significant  effect  on  the  turbulence  level  inside  the  sensor.  This  is  believed 
to  be  due  to  separation  of  the  flow,  see  Figure  24.  The  small  radius  of  curvature 
model  generates  turbulence  inside  the  sensor  from  the  separated  flow  region.  It 
is  significant  that  the  separated  flow  only  occurred  at  the  bottom  half.  This  is 
believed  due  to  the  fact  that  distance  between  the  gap  and  the  top  of  the  instrument 
is  two  inches  whereas  the  mounting  system  resulted  in  much  longer  cylinder  be- 
neath the  corona  gap.  This  results  in  a thicker  boundary  layer  on  the  bottom  and 
produces  the  separated  flow.  The  corona  instrument  itself  effectively  integrates 
the  turbulence  along  the  corona  path  to  provide  a measure  of  the  average  turbulence 
in  the  corona  sensor.  The  averaging  reduces  the  effect  the  separated  flow  can  have. 
The  observed  power  spectral  density  data.  Figure  21  demonstrates  that  separation 
induced  turbulence  does  not  significantly  alter  the  spectral  shape  at  turbulence 
levels  greater  than  one  percent.  Nevertheless,  care  should  be  taken  to  increase 
the  leading  edge  radius  in  order  to  minimize  the  spectral  contamination  from  locally 
generated  turbulence. 

The  corona  instrument  was  calibrated  for  velocity  by  comparing  known  wind 
tunnel  velocity  with  sensor  voltage  output.  The  results  are  shown  in  Figure  25. 

The  sensor  was  tilted  in  and  away  from  the  wind  vector  up  to  30  degrees.  At  no 
time  during  the  tilting  was  the  measured  velocity  inside  the  sensor  observed  to 
change. 
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Figure  25.  Wind  Tunnel  Calibration  of  Velocity  Including  Effects 
of  Angle  of  Attack 


3.  I INITIAL  FLIGHT  TEST 


Two  CADT  sensors  were  constrained  to  a 10  ni  vertical  separation.  This 
allowed  measurement  of  wind  shears  and  cross-correlation  coefficients  of  the 
vertical  wind  profile.  Once  launched,  the  electronics,  telemetry,  and  sensors  are 
reeled  down  200  m below  the  balloon  load  bar,  thus  eliminating  interference  from 
the  balloon  itself.  The  two  sensors,  protected  in  an  aluminum  gondola,  are  ex- 
tended outward  and  atmospheric  data  are  gathered  from  40,  000  to  100,  000  ft.  The 
payload  remains  aloft  for  10  to  12  hr  and  the  balloon  is  positioned  to  altitudes  of 
maximum  interest.  Data  are  reduced  and  analyzed  by  a real  time  computer  to 
yield  values  for  mean  differential  winds,  wind  shears,  and  turbulent  intensities. 
These  values  are  correlated  with  space  and  time  to  allow  construction  of  the  ver- 
tical wind  profile  and  turbulent  layering.  We  are  able  to  watch  the  historical 
development  of  clear  air  turbulence  which  contributes  to  models  for  stratospheric 
transport. 

Since  the  ion  beam  originates  from  a random  spot  on  the  needle  and  slowly 
travels  about  the  needle,  the  zero-wind  position  of  the  ion  beam  centerline  as  it 
falls  on  the  target  plate  must  be  determined.  This  is  done  by  periodically  closing 
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a shutter  over  the  airstream  gap.  The  sensor  automatically  measures  the  zero- 
wind  offset  currents  and  corrects  it  to  zero. 

Upon  completion  of  the  experiment,  the  sensors  are  reeled  back  into  their 
protective  shrouds.  The  payload  is  cut  from  the  loadbar  and  descends  on  its  own 
parachute  to  a site  where  a recovery  crew  is  standing  by. 

The  initial  flight  was  conducted  from  Holloman  AFB  on  28  July  1976  at  4 AM. 
The  corona  instrument  was  turned  on  at  15.  1 km  with  measurement  of  stratospheric 
wind  shear  and  turbulence  conducted  up  to  15.8  km.  The  experiment  was  termi- 
nated after  10  min.  of  data  due  to  loss  of  power  resulting  from  damage  to  the  instru- 
mentation package  and  one  sensor  during  launching  of  the  balloon  payload.  The  alti- 
tude excursion  measured  in  this  brief  flight  is  shown  in  Figure  26.  The  measured 
pressure,  temperature,  and  velocity  sensed  by  the  corona  instrument  is  shown  in 
Figure  27.  The  velocity  measured  represents  the  velocity  difference  between  the 
balloon  and  the  corona  anemometer.  The  anemometer  was  reeled  down  296  m from 
the  balloon  to  avoid  balloon  wake  effects.  It  is  assumed  that  the  172  ft  diameter 
balloon  drifts  with  the  wind. 

The  temperature  minimum  occurred  at  15.45  km.  Thus,  this  first  flight  yields 
a measurement  above  and  below  the  tropopause.  The  turbulent  intensity  of  this 
region  as  determined  from  the  variance  of  the  velocity  is  shown  in  Figure  28.  A 
large  turbulent  layer  is  located  just  below  the  turpopause.  Spectral  power  density 
distributions  were  computed  from  the  measured  velocity  and  are  shown  in  Figure 
29.  The  data  show  the  flight  instrument  has  a noise  power  spectrum  distribution 
with  a -5/3  slope  with  frequency.  The  first  spectra  shown  for  time  0 exhibits  the 
same  -5/3  slope.  As  the  balloon  rises  and  passes  through  the  large  turbulent, 
layer,  the  sequence  of  power  density  spectral  indicates  a rise  and  then  decay  in 
spectral  energy  near  200  Hz.  An  interesting  feature  is  the  -6  slope  of  the  spectral 
energy  that  exists  when  the  corona  is  in  the  turbulent  layer.  An  altitude  profile  of 
turbulent  dissipation  rate  is  given  in  Figure  30.  This  curve  is  very  similar  to  the 
velocity  variance,  as  should  be  expected. 

The  results  contained  in  this  report  are  given  to  illustrate  the  operational 
capability  of  the  corona  anemometer.  A thorough  evaluation  of  the  measurements, 
including  analysis  of  the  rate  of  turbulent  dissipation,  evaluating  of  turbulent  scale 
lengths  and  characterization  of  the  turbulent  layers  will  be  presented  in  another 
technical  report. 
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Figure  27.  Observed  Atmospheric  Pressure  and  Temperature  and  the 
Velocity  Measured  by  the  Corona  Anemometer 
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Figure  30.  Altitude  Profile  of  Turbulent  Dissipate  Rate 


4.  CONCLUSIONS 


A new  wind  and  turbulence  sensor  capable  of  operating  in  the  stratosphere  has 
been  designed,  calibrated  and  flight  tested.  The  sensor,  utilizing  the  principle  of 
a corona  ion  beam  deflected  by  the  airstream,  is  shown  to  function  properly  and 
measure  turbulent  wind  velocities  out  to  a frequency  of  2000  Hz.  The  corona  sensor 
data  can  be  analyzed  to  directly  determine  the  turbulent  energy  contained  in  the  thin 
turbulence  present  in  the  stratosphere.  Two  sensors,  separated  some  10  m can  be 
flown  to  directly  measure  the  vertical  gradient  of  the  horizontal  wind  velocity  and 
the  temperature  gradient.  This  permits  a direct  evaluation  of  the  Richardson  num- 
ber, a parameter  related  to  the  occurrence  of  turbulent  layers  in  the  stratosphere. 
Future  flights  are  planned  to  obtain  altitude  profiles  of  the  turbulence  layers  and 
the  time  history  of  given  turbulent  layers.  The  latter  can  be  obtained  by  position- 
ing the  balloon  to  float  in  a turbulent  layer. 

As  of  28  April  1978,  six  individual  flights  have  taken  place  from  Holloman  AFB, 
NM.  Refined  methods  of  calibration  and  analysis  have  been  developed.  Also  a new 
ion  tlme-of-flight  instrument  has  been  developed  which  can  make  in-situ  measure- 
ments of  the  DC  winds.  These  developments  will  appear  in  another  report. 
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